Ferromagnetism and orbital order in a topological ferroelectric 
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We explore via density functional calculations the magnetic doping of a topological ferroelectric 
as an unconventional route to multiferroicity. Vanadium doping of the layered perovskite La2Ti207 
largely preserves electric polarization and produces robust ferromagnetic order, hence proper multi- 
ferroicity. The marked tendency of dopants to cluster into chains results in an insulating character 
at generic doping. Ferromagnetism stems from the symmetry breaking of the multi-orbital V system 
via an unusual "antiferro" -orbital order, and from the host's low-symmetry layered structure. 



PACS numbers: 75.85. +t,77.84.-s, 75. 30.-m,71.15.Mb 

Multiferroic materials, where ferroelectricity and mag- 
netism coexist, have been studied intensively in recent 
years [1] due to their fascinating physical properties 
and their potential for application (e.g. multi-state or 
electrically switchable magnetic memories) [2]. Despite 
the effort, single-phase bulk multiferroics at room tem- 
perature are still unavailable; most are antiferromag- 
netic (AF) and not, as desired, ferromagnetic (FM); 
and, generally, they do not meet [3] the integration 
and functional requirements of device applications [H [5] . 
A parallel line of research [31 H] on multi-phase sys- 
tems (interfaces, heterostructures, nanostructures) re- 
vealed unconventional ferroelectricity mechanisms in ar- 
tificially layered perovskites (e.g. SrTiOa/PbTiOs super- 
lattices [7|). Recently, it was suggested [5J [§] that mag- 
netism and ferroelectricity will be more cooperative in 
lower-symmetry, naturally nano-structured, layered sys- 
tems. In one such system, the layered perovskite ox- 
ide La2Ti2C>7 (LTO), ferroelectricity is not displacive 
but rather results from dipoles generated by antiferrodis- 
tortive oxygen-octahedra rotations, which fail to cancel 
out due to the layered structure, resulting in a net macro- 
scopic polarization [5]. Because ferroelectricity follows 
directly from its layered structure, LTO was labeled [9] a 
topological ferroelectric. LTO has a Tc of 1770 K, a siz- 
able polarization P c —h /xC/cm 2 jTUj along the c crystal 
axis (schematic in Fig. 2), and non-critical permittivity 
e~50 which would be an asset for integration in ferro- 
electric random-access memory elements jllj. 

Here we explore theoretically a route to multiferroic- 
ity based on doping LTO with magnetically active ions 
substituting for Ti (nominally a 4+ 3d ion). By an ex- 
tensive screening of transition-metal dopants, we found 
vanadium (V, nominally 4+ and 3d 1 in LTO) to yield 
robust ferromagnetic order while preserving polarization 
and a small gap, thus giving rise to proper multiferroicity. 
We find that V align in chains along the a axis, orthog- 
onal to the polarization c axis (with magnetic moments 
orthogonal to the chains themselves), resulting in insu- 
lating character at generic doping. Intriguingly, ferro- 
magnetic coupling is linked to the layered structure. The 



latter "pre-wires" a two-orbital-per-site system, whose 
degeneracy is broken in favor of ferromagnetic superex- 
change by a peculiar orbital ordering, as foreshadowed 
in earlier theoretical work |12j . We preliminarily investi- 
gate the magnetoelectric response to (a subset of) ionic 
vibrations, finding it to be almost absent. 

We perform ab initio calculations within spin-density 
functional theory (DFT) in the generalized gradient ap- 
proximation (GGA) [T3] as implemented in the VASP 
code [14], with projector-augmented waves [15] for La 
5spd6s, Ti and V 3pdAs, and O 2sp states, plane-wave 
cutoff 330 eV, and a 4x2x5 grid in the primitive cell's 
Brillouin zone, folded for larger supercells. Magnetic 
calculations for magnetoelectricity and anisotropy are 
non-collinear and include spin-orbit coupling (see below) . 
GGA seems a reliable approach for LTO, which is a band 
insulator with O p valence and Ti d conduction [5] ; even 
the V electrons (see below) sit into conduction-derived 
states, whose the localization is considerable but not ex- 
treme. 

LTO is the n—A member of the A n B n 03 rl +2 family, 
and consists of four perovskite-like units (distorted Ti06 
octahedra) stacked along the [Oil] cartesian direction 
and separated periodically by the insertion of an extra 
O layer [TB]. We optimize cell and internal parameters 
(force threshold 20 meV/A) in the experimentally-stable 
monoclinic P2\ structure. Our calculated lattice param- 
eters a=7.80 A, 6=13.22 A, c=5.58 A, and 7=98.52° are 
in excellent agreement with experiment |17j . The ener- 
getics and electronic structure of V-doped LTO is ob- 
tained collating data from calculations on 44-atom and 
88-atom cells at 12.5% average V concentration. Our 
conclusions do not really depend critically on this fig- 
ure, which in essence only amounts to approximate the 
isolated-impurity limit with a single V at no less than 
about 8 A from its periodic image. The isolated V's spin- 
polarized electron yields a magnetic moment of 1.0 /xb- 
Its spin-resolved and atom-projected density of states 
(DOS) shown in Fig[l] is typical of all FM configura- 
tions of V discussed below. The valence band is mostly 
of O p character, and the conduction band is mostly of 
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O-hybridized Ti d character. Near the Fermi level, V 
induces two spin-majority peaks, one carrying the extra 
electron and the other empty, across a gap of 0.3 eV. Ti 
states start 0.5 eV, and minority V d states at ~0.85 eV 
above the edge of the occupied V band. 
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FIG. 1. Total density of states (DOS) of LTO:V (gray shade) 
and its site projections for O (dashed red line), Ti (solid blue 
line), and V (green shade). 

To estimate the magnetic coupling between V's, we 
sampled a variety of structural (V pairs on the relevant 
non symmetry-equivalent cation B sites) and magnetic 
(collinear FM and AF) configurations. Structurally, we 
find a tendency to chain ordering of V along the a di- 
rection (see Fig{2]). Chains are favored by 100 meV per 
V pair over the first competing configuration, and 180 
meV per V pair referred to isolated impurities. Magnet- 
ically, FM order is very stable within the V chains, with 
FM couplings J a ~80 meV and 20 meV along the a direc- 
tion, respectively, for chains in proximity of the extra O 
layer or inside the perovskite block. Both types of chain 
(or fragment thereof, see below) will appear in similar 
concentrations as their energies differ by only about 5 
meV/dopant. The interchain coupling in the b,c plane is 
also FM and small, J&, c ~3 meV at 12% V. 

For a ballpark estimate of the critical temperature we 
use the Ising model with axial anisotropy [IB], whose 
critical temperature Tc is determined by the ratio 
R=Jb.c/Ja of in-plane to on-axis coupling. The model 
mirrors closely our system of weakly coupled chains with 
strong internal coupling along one axis. We estimate 
the FM Tc for both our two nearly-degenerate ground 
states with distinct J a 's, since magnetic couplings are 
not additive. For J a ~80 meV, i?=0.04 and correspond- 
ingly T c ~0.93J a =860 K; for J a ~20 meV, i?=0.15 and 
Tc~1.55 J a =360 K, i.e. Tc is around room temperature 
in the worst-case estimate. 

We now check that the material remains ferroelectric 
upon doping. The polarization calculated via the Berry 



phase technique [H] for LTO is P°=5.2 /zC/cm 2 along 
the c axis, similarly to experiment and previous calcula- 
tions [TO] [20]. In V-doped LTO, the distortions gener- 
ating ferroelectricity are fully conserved, aside from mi- 
nor quantitative differences. We find a slightly enhanced 
and reduced P c for the energetically-favored chain con- 
figurations, respectively, inside the blocks and near the 
additional O layers. Since polarization is extensive, we 
average the values for these nearly-degenerate configura- 
tions obtaining P c =4.7 /zC/cm 2 , i.e. only 10% lower than 
undoped LTO at 12.5% V concentration. From P values 
at 12% and 25% V, we infer that P c [x]~P° '(1-x) with x 
% V. Interestingly, V doping produces polarization com- 
ponents along a and b, 4.5 fiC/cm 2 and 0.7 fiC/m 2 re- 
spectively at 12% V, comparable with that along c. We 
will address the consequences of this elsewhere. 

So far we saw that V-doped LTO is insulating, polar- 
ized, and ferromagnetic, and therefore properly multifer- 
roic. To show that this applies at generic doping, we 
evaluate energetics, magnetism, and insulating character 
of V chain fragments in LTO comparing three configura- 
tions with 4 cations lined up along a: i) -V-V-V-V-, ii) 
-V-V-Ti-Ti-, and in) -V-Ti-V-Ti-. The full chain i) is 
bound by 100 meV compared to two separated V dimers 
(chain fragments) as in ii). In turn, the V dimer frag- 
ment ii) is 130 meV lower in energy than the "broken" 
dimer Hi). Energy thus drives V to cluster in homoge- 
neous chains; these will be in fact finite chain fragments 
accounting for configurational entropy (we will elaborate 
on this elsewhere.) All chain and chain-fragment config- 
urations are insulating, and i) and ii) have similar FM 
coupling along a as measured by the FM-AF energy dif- 
ference; these are clear indications that the system will 
be insulating and FM at a generic concentration due to 
the V tendency to chain clustering. 

We now analyze the strong ferromagnetic coupling. V- 
doped LTO is at first glance a degenerate multi-orbital 
system. However, this degeneracy is resolved so as to 
favor FM order, and the mechanism is, to a significant 
extent, built-in into the host structure. First, only two 
orbitals per V site are involved: the V 3d DOS of a full 
V chain in Figj3^, shows d xy , d xz states lower in energy 
and involved in V-V bonding and magnetism, while d yz 
follows the e g states to higher energy. Indeed, in the 2' 
point group of LTO the octahedral-symmetry ti g and e g 
representations split into 1-fold a\ and b\. d yz and e g 
functions belongs to a%, whereas d xy and d xz belong to 
b\ . Physically, one notices that the periodic alignment of 
octahcdra along their axes is interrupted by the layered 
structure, except along the a axis on which d xy and d xz 
overlap significantly. Second, the bonding-antibonding 
combination (Fig[3]) of d xy and d xz along the chain gen- 
erates two narrow bands separated by a small gap (see 
e.g. [H] for a similar case). Third, because of the large 
energy gain due to Hund on-site exchange, FM order will 
clearly be favored over AF (whose much smaller gain is of 
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FIG. 2. On the right: Isosurface plot of the spin density for 
the LTO:V in V-chain configuration; on the left: 2D cuts of 
the spin density along the two A and B planes, normal to the 
a axis, in proximity of the VI and V2 atoms, respectively. 



order the superexchange energy), provided that, fourth, 
the two-orbital bonding combination is non-degenerate, 
and hence different on each V site, so as to host two spin- 
parallel electrons. Indeed, Fig(3^t clearly shows the dif- 
ferent weight of the two orbitals in the bonding state on 
a specific site; this alternation is accompanied by mod- 
ulated bond lengths and by differences in V magnetic 
moments, as shown below. The resulting FM spin order 
and "antiferro" -orbital order for the orbital filling agree 
with early predictions [12] as well as with Goodenough- 
Kanamori rules [35]. The chain-like ordering of V, we 
note, is actually relevant even in some binary vanadates, 
as discussed e.g. in [2"3"] . 

The computed spin density shown in Fig[2]further con- 
firms the above argument. It is localized on V atoms and 
has the fingerprint shape of a i2g-like combination. The 
2D cuts of the spin density along planes (A and B in the 
figure) normal to the a axis and close to the V atoms, 
clearly show "antiferro" -orbital order. The spin density 
lobes on the two V sites are distinctly slanted towards 
alternating O atoms in the upper layer, one roughly nor- 
mal to the other. The V chain has resolved the degener- 
acy by alternating d xy -d xz orbital combinations a and /3, 
respectively with 65/35% and 54/46% d xy /d xz content. 
As hinted-to above, concomitantly with spin and orbital 
ordering, the V-chain is also characterized by dimerized 
bond lengths and moments (see Fig (2]): pairs of short 
(1.95 and 1.97 A) and long (2.00 and 2.02 A) O-V-0 
bonds alternate along the chain direction, and the mag- 
netic moments of short-bond and long-bond V are also 
alternating (0.75 and 1.03 hb)- As seen in FigJ5J state /? 
is centered on the long-bond large-moment V2 atom and 
a on Vi. 
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FIG. 3. (a) Projected orbital resolved DOS relative to a V 
atom in LTO:V with 2V/cell, the V-chain configuration, (b) 
Schematic energy levels diagram for interacting V-V pair in 
LTO:V. 



To further validate the link between V dimerization 
and orbital ordering, we compare V-doped LTO in P2\ 
symmetry with the parent high-symmetry structure with 
space group P2\/m. This structure has no antiferro- 
distortive rotations and only one inequivalent V site 
along each a chain, with V-V distances 3.91 A and 180° 
O-V-0 angle. Thus, it cannot remove the orbital degen- 
eracy, and the bonding-antibonding splitting does not 
materialize. The system ends up being metallic and fer- 
romagnetic, with identical spin-density on all V atoms. 
The structural dimerization is evidently an integral part 
of the explanation of the insulating nature and the spin 
and orbital ordering in V-doped LTO, as its absence 
causes orbital degeneracy and metallicity. In this very 
light, and in keeping with our previous symmetry con- 
siderations, we note that the O-Ti-0 bond lengths along 
a are already modulated in LTO (1.97, 2.03, 1.98, and 
2.02 A). The structural inequivalencies along the a axis 
are thus at least partly pre-wired in the host structure's 
low symmetry, supporting the general idea of a search 
for multiferroicity opportunities in unconventional, low 
symmetry materials. 

Since multiferroicity and magnetoelectricity (the cou- 
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pling of, e.g., electric field to magnetization) are often 
associated, we study preliminarily the lattice compo- 
nent of the magnetoelectric tensor, expressed by Ihiguez 
[24j as the product of polarization and magnetization 
changes in response to polar-mode distortion divided by 
the mode force constant. We need magnetic non-collincar 
spin-orbit calculations with full spin relaxation and a 
stricter convergence limit than usual (<0.01 fieV). Wc 
consider small displacements away from the ferroelectric 
P2i structure along the composite IR mode (mainly rota- 
tions around a) connecting it to the paraelectric P2\jm 
structure. To numerical accuracy, all magnetization com- 
ponents are insensitive to this composite displacement. 
A rationale is provided by magnetic anisotropy calcula- 
tions, whereby we constrain the spin direction to map 
out the energy vs spin orientation. The V spins prefer- 
entially point normal to the chain direction-that is, a is 
the hard axis, with a substantial anisotropy energy 180 
/xeV/ spin with respect to the b,c plane. The b,c plane, in 
turn, hosts many spin-configurational minima within less 
than 5 /xeV/spin, i.e. we all directions in the 6,c plane are 
equally easy: thus, spins "freewheel" around the a axis ir- 
respective of the rotations (we dub this the "ball-bearing" 
effect) and are insensitive to octahedra rotations around 
a. Interestingly, V spins are always approximately par- 
allel, and they point near the [Oil] direction in the b,c 
plane in the lowest-energy configuration. In conclusion, 
for the composite ferro-para IR mode the lattice mag- 
netoelectric tensor is zero. This may not be the whole 
story, however, as we did not consider all possible IR 
modes; some of those we neglected may couple to the 
a and b polarization components appearing in V-doped 
LTO. Also, we ignored the electronic component, which 
might still be appreciable. 

Finally, we investigate the effect of oxygen vacancies 
(which are donor defects) on magnetism, motivated by 
our finding for Sc-doped LTO that the vacancy forma- 
tion energy drops to almost zero, destroying ferromag- 
netism. In equilibrium oxygen-rich conditions (chemical 
potential /io=A 1 o 2 /2), the vacancy concentration is neg- 
ligible at typical growth temperatures, because the low- 
est (as function of configuration) formation energy of the 
vacancy in the neutral state in LTO is 5.4 eV. In V- 
doped LTO, the concentration is still negligible, as the 
formation energy is 4.6 eV. Consistently with the lower 
formation energy, the vacancy binds to V pairs within 
the a chain, turning the magnetic ordering from FM to 
AF. This suggests the possibility to spatially structure 
the magnetic state by controlling the local O content in 
layer-by-layer growth techniques. 

In summary, the V-doped layered ferroelectric LTO 
shows several unique properties of basic as well as ap- 
plicative (e.g. data storage) interest. It is a ferromag- 
netic multiferroic with high temperature ferroelectricity 
and ferromagnetism at any doping, due to V clustering in 
longish, homogeneous chains. Polarization is significant, 



and decreases mildly as function of V concentration (at 
least at small doping). Magnetoelectric coupling seems 
negligible at the present approximation level. Control of 
native defects could turn FM order into AF, enabling the 
engineering of resistivity and magnetic phase changes via 
defect injection and migration. 
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